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ABSTRACT 


The perturbations of alternating geomagnetic fields 
of global dimensions by two-dimensional local structures is 
studied in a cylindrical co-ordinate system. A general 
numerical finite difference technique is applied to solve 
Maxwell's equations,at,each.point.of a polar, mesh »with 
variable grid spacings. 

The proper size and spacing of the mesh used for 
different conductivity configurations is discussed with 
respect to the Gauss-Seidel iterative method which is used 
to solve the finite difference equations. 

Two symmetric models, one a ridge-like structure 
andrthe \otherwaedyke, areistudied for |\differentscoscillation 
periods of the source. A series of sloping contacts with 
different slopes is investigated as well. In the H-polariza- 
tion case the magnetic field and in the E-polarization the 
electric tteld ais given over the non-uniform region at 
different epochs of the oscillation period. Also, the 
electric and magnetic field components as well as the 
apparent resistivity and the) ratio of the vertical and 
horizontal magnetic components are calculated along the 


surface. These profiles are discussed in detail. 
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CHAPTER 1 INTRODUCTION 


Pol as torical Background 


Many studies have been made to determine the 
pertiuirbati ons tof #al ternatins: electromagnetic Fields shy 
dtfterent geological conf feurati ons ; Some authors have 
approached tthe problem analytically (d'Erceville and Kunetz, 
1962; Weaver, 1963; Blake, 1970; Weaver and Thomson, 1972; 
Geyer, 1972) while others have made analogue model measure- 
ments (Neves, 1957; Dosso, 1966). Several studies have 
applied numerical methods to obtain solutions for Maxwell's 
equations. Some of the most often used of these latter 
methods are the finite difference technique (Neves, 1957; 
Patrrvel Cand “Bosttwens L969) ones "and *P niceln «19 70). 19 7 bax OOD = 
Jones fand *Pascoe ,lO7/ 1S sPastcoe stand Jones, (b9r7-2)if tehe 
transmission “liane analogy (Madden “and Thompson, 91965); 
Wright, 1969) and the finite element method (Coggon, 1971; 
Reddy. sand. Rankin ¥ 91917/3)% 

Since ‘many “geological structures fare ofitwo- 
dimensional nature, such an approach to the geomagnetic 
induction problem has often been taken, and in most of 
these works the Cartesian co-ordinate system has been used. 
In some cases, especially when sloping contact models are 


studied, it seems more reasonable to use a cylindrical 
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co-ordinate system. Geyer (1972) considered analytically a 
stoping contact in cylindrical co-ordinates yibitviound wit 
necessary to assume that the tangential component of the 
Magnetiewiieldsin the E-=polarization case is sconstanthabsh 
along the surface. However, as discussed by Jones and 
Renee (1970) this *approximation is not adequatesto gdeseuibe 
Ehe@tgtetnaturevor the fields. 

The finite element method using triangular meshes 
is one approximation to the sloping contact model. Reddy 
and Rankin (1973) gave some general conclusions concerning 
the electric and magnetic field values along the surface for 
auconvtact which|is initially sloping but becomesthorizontal, 
However, in their paper very little detail of the mesh 
dimensions®on number of -gridpoints is given. Lteis »thereftore 
appropriate to investigate some of these conclusions further, 
and#ingperticubar with» reference sto a purely sloping? contact. 

It is reasonable to approach the problem by applying 
theptinitegdiftfierence=method in agcylimdrical co-ordinate 
system. atlhis is suitable for a purely sloping contactgand 


can also be used for other models of interest. 


ii2) General Outline of ‘the Study 


The problem is to determine the electric and 
Magnetic taields associated with currents induced in the 


Earth by given geomagnetic sources. When limited areas are 
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concerned, ThhesEarthecanobél treatedhasaaesemi=infinite 
composite conductor which may have horizontal layers and 
lateral inhomogeneities and has a plane surface. This 

local problem may be considered as a perturbation problem 
(Beicegtdoaeyid The thtuced&field.is ere ned by the 

focal redistributionyofea currentesystemhof hargei dimensions? 
tihte globaisindtedngl fiekdmisauniformrandt}has*thensame 
intensity and form all over the region considered. In the 
present study the source field shows sinusoidal time 
Variation. 

The work of this thesis is concerned with the 
numerical solution. of the perturbation problem in cylindrical 
co-ordinates where the conductivity configuration does not 
change in the z direction. This implies that the problem is 
of two-dimensional nature and a polar mesh is used. Maxwell's 
equations are solved at each gridpoint applying the finite 
difference technique. The method is general so that the 
sloping contact as well as other conductivity configurations 
may be studied. in particular a buried cidge—-like upwelling 
(Figure 8) and a dyke-like structure (Figure 9) as well as 
five different sloping contacts (Figure 10) are examined 
here. The electric and magnetic field distributions within 
and outside the conducting region are studied as well as 
the nature of the current flow associated with the inhomo- 


geneities for two polarization cases. 
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The surface profiles of the electric and magnetic 
fields are of particular interest since the field components 
induced by ionospheric and magnetospheric sources may be 
measured above natural geological configurations by the 
magnetotelluric and other methods. Large peotoeical structures 
which exist naturally can be recognized by comparing observed 


data with computed electromagnetic perturbation models. 
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CHAPTER 2 MATHEMATICAL FORMULATION OF THE ELECTROMAGNETIC 
PROBLEM IN A CYLINDRICAL CO-ORDINATE SYSTEM 


Zo spasic Bquati ons 


The co-ordinate system is shown in Figure l. 
Maxwell's equations are solved in this system where the con- 
ductivity configuration and all field values are constant in 
the z direction. Since displacement currents may be ignored 
for the frequencies of the geomagnetic variations considered 


(Lahiri and Price, 1939), Maxwell's equations become 
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i} 


curl H koE C22) 


where it is assumed that the time factor for all field 
quantities is, exp (iwet). In electromagnetic units wu = 1 and 
k = 41, while in MKS units u = 4m x NG cupeardakew lat 

As E and H are independent of z, the six equations 
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Figure 1 The cylindrical co-ordinate system. 


Figure 2 An arbitrary part of the mesh surrounding 
gridpoint ‘O°; 


Figure 3 The Location of gridpoints near the surface. 
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Equations (2.dce), (22a) and (2.25) involve only 


Ho E. and ay and represent the H-polarization case, while 


Gquations. (2.1la)5 (2. ib) and (2e2¢) anvolve E ee and Hy 


and represent the E-polarization case. These two sets of 


equations can be solved independently. The elimination of 


ES. and . from equation (2.1c) gives 
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for the H-polarization case, while eliminating He and H 
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from equation (2.2c) gives 
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for the E-polarization case. 
If F is used to represent H, or E. depending on the 


case being considered, the same equation may be used for 


equations (2.3),and (2.4) 
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where no = ukwo , 
Since the field quantity F = f + ig is a complex 


number, equation (2.5) may be written as two equations 


Vi=-n 8g Cray) 


Ue eusen ce (2.5b) 


Equations (2.54) and (2:56) must be solved 
simultaneously in all regions being considered. In addition, 


the boundary conditions have to be satisfied. 


2.2 The Numerical Method 


In the numerical method the conductivity distribu- 
tion is given over a two-dimensional polar mesh, whose 
radius is great enough to satisfy the external boundary 
Condittons:y WAN arbitrary part*or this*¢eriad “is shown in 
Figure 2% 

In the present work equations (2.5a) and (2.5b) 
are replaced by finite difference equations. Any point, 


except those along the external boundary, can be chosen as 
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Taylor's theorem is used to obtain the correspond- 


ing finite difference equations. 
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There is a similar set of equations for the g values. 
If higher order terms are neglected and finite 
differences. are substituted for. the first derivatives, 


equations (2.6a) and (02.66) can be written: 
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the £irst derivatives=ot functions £ and <¢ may be 
approximated by finite differences as above if they vary 
nearly linearly in the small distances concerned and the 
adjacent grid spacings do not differ from each other very much. 

In the present work the Gauss-Seidel iteration 
method is employed wherein the most recent values of f and g 
for the four nearest points surrounding point '0Q" are used in 
each step of the iteration to determine the subsequent values 
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2.3 Internal Boundary Conditions 


2,51 General 


Equations (2.8a) and (2.8b) must be satisfied at 
point '0' for each of the regions surrounding it. This means 
that both (2.8a) and (2.8b) represent four equations each 
since there are four regions of different conductivity around 


point '0!. ,These equations in which the firsce subscript 
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refers to the conductive region while the second subscript 


indicates the grid point are written 
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The underlined values are so-called ‘fictitious values' 
(Jones and Pascoe, 1971) and must be determined from the 


boundary conditions. 
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existing values. 


The boundary conditions are 
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Theghictitioussvalues of Bes = Ete + 185% may be 
€kininated gingsequatieons)(2.9.a)eand «(2 <9b). enSinee tye = f. 
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in thest=polardizatione caséeecequatd onsei(24i1da) and 
(2.11b) are solved simultaneously at each interior point of 
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By “substitution of “the underlined fictitious values 
from the boundary conditions into equations .(2 294) and (2.96) 


after adding we obtain a pair of equations 
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where A, b, Bi» B,> B and By are as before. 
In the E-polarization case the values fo and Bo 


at each interior point are determined from equations (2.13a) 


and, (2.15b) usdne Cramer's Rule. 


2.4 Boundary Conditions at the Surface 
2.4, H-polapization 


In the H-polarization case, it is seen from 
editions 1G2. 2a) sand “(2 -2b)) that H is constant everywhere 
OuUreiae the conductinge repion. ssince sale components soln 
are continuous across any interface, this implies that all 
along the surface the value of H. is constant and equal to 
its value in the non-conducting upper half-space. The surface 
of the conducting region then becomes one of the external 
boundaries for the H-polarization case. Also, the magnetic 
field &@s normalized such’ that 1 rd. (1.04 9050) ts the valee 


along this boundary. 
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In the E-polarization case there is no correspondingly 
Simple boundary condition for E Price (1950) pointed ,out 
that the problem is completely determined if the two-dimen- 
sional half-space is considered as the limit of a three- 
dimensional spherical conductor. In” this case the ratio of 
the tangential components of the inducing and induced fields 
ie independent of the conductivity as the, vadiuse? tends (Ge 
infinity. Since the inducing field is assumed to be of the 
same intensity and form over the entire composite conductor, 
He is the same at the surface approaching the external 
boundaries (Jones and Price, 1970). At these extremities 
of the mesh ma is near zero if the boundary is far enough 
from any discontinuities inside the conductor near tne 
Burtace, in this region the field configuration may be 


regarded as that of a uniform conducting half-space and the 


Sorutton of equation C2,4) as 
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where E is the value of E. at the surface. From equation 
Oo 


(2.la) we obtain the surface value of He 
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As r tends to infinity the radial component equals the total 


value of the magnetic field along the surface and 
H =H =-H C2ioay 


wheres thie, indices, Rand. L refer to, the values at the right 
and left hand sides of the model respectively. From equation 


(2.15) we obtain the boundary condition for ae 
‘ (2415 b) 


The values Np and ny are determined by the conductivities at 
the extremities of the model just inside the conductor. The 
electromagnetic field is normalized such that EOL SA Gl a0 00) 
is taken for the computations. 

In the E-polarization case the field values must 
be calculated for the non-conducting half-space as well as 
im,che sinterior of the conductor’ Equation (2.2c) indicates 
that HL is constant at the extremities of the model near the 
surface in the non-conducting region. E. here can be deter- 


timed irom ~equations (2415) .and (21a), since all oehe wtangen— 


tial components of E are continuous across any interface. 
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CHAPTER 3 GENERAL DISCUSSION OF THE COMPUTATIONS 


3,14; Initial.Values.for’the Iteration 


in tie H-polarizgation case Gquations (2.1lia ana Dp) 
and in the E-polarization case equations (2.13a and BD) are 
simultaneously solved for f and g applying the Gauss-Seidel 
iteration method (Smith, I969)°" The way in which the initial 
values of the unknown quantities are set over the mesh in- 
Fluences the number of iterations necessary to obtain the 
final values. If the initial values exhibit the expected 
variations in the field configuration, the process requires 
fewer iterations. 

In the H-polarization case the initial field values 
are set using the solution of equation (2.03) in 4 Uniform 
conducting half-space 


Ue ge 
Z 


NG sim oy i ‘Cena, 


He a H exp (- 


Ho is the normalized surface value and n refers to the 
conducting region where the actual gridpoint is located. 
Ti the point ts on thé boundary of regions of different 
conductivities, then the average value is used for n . 

There is also another way to set the initial values. 
tt is possible to use the appropriate fj values in equation 
(3.1) only at great distances from the discontinuities where 


the field may be regarded as uniform and to apply a linear 
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interpolation for n between the external boundaries in the 
rest of the conducting region. A vertical contact model was 
studied by setting the initial values in both ways. The 
number of iterations using the latter set of initial values 
significantly exceeded the number of iterations using the 
actual values of 1) in/e@ach region. It is apparent that the 
sudden changes in the field values near the vertical discon- 
tinuity are smoothed by the linear approximation in n and it 
takes longer to approach the real field configuration. 

In the E-polarization case the initial field values 
are determined from equation (2.14) for the conducting halt= 
space and from equation (2.16b) for the non-conducting region. 
ivequation (2 .L6b). the. quantlties. jn iand E represent inter- 
polated values’, A linear change across the upper half of the 


grid is assumed for (n n1? as well as for (EoR a ae 


R oL 


Similarly, as in the H-polarization case there are 
different ways to set the initial values in the bottom half- 
space by choosing n and E in equation (2.14). Both of them 
may be determined using a linear interpolation as in the 
non-conducting region. Also in the E-polarization a vertical 
contact model was examined by setting the initial f and g 
values in different ways. The number of the iterations 
required to approach the final values is the smallest when 


the actual n values are used in each grid point and there is 
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a sudden change in E_ at the discontinuity. The applied E 
O ) 


values corresponding to equation (2.15b) are 


Set for vices 4 (3328) 
Up 
Eo oe ir for NV = Np (3525) 
if} 
ee ee 
UR 
E a aarang ewe date thepdiscontinuity . Ose Ley 


The iteration procedure takes considerably longer when a 
linear approximation is used either for n or Eo or for both 
of them. 

In other models the initial EG values are determined 
fromeequatione (3.20, bo: and ¢).) ‘Equation Mo.2c) mey, Desused 


at every gridpoint where n does not equal is nor Np 


Sele Loe Grad 


The meshes used for the two different polarization 
cases are illustrated in Figure 4. The grid spacings are 
variable so that the values Ag, (4 = l,m) and Sees = l,n) 
may vary. Also, each region of the mesh can have material 
of different conductivity. To satisfy the boundary conditions 
there must not be conductivity changes near the surface in 
the outer Ar regions. Furthermore, in the E-polarization 


case o=0 for 180° < od < 360° 
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The distancetirom thesorigin étotthe 4external 
n 


boundaries Gh s6 Some Ar,) is determined for each model 
J=1 
considering the skin depths and the conductivity configuration. 


The «skin depth is calculated foreach (material sin ithe model 
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- leas 5 (353) 


The field values along the arc of maximum radius do 
notiichange during the iteration procedure is Therefore f ithe 
initial values of the iteration must be good approximations 
of ithe ‘real 'fveld values there: THis implies 'that athe 
external boundary must be far enough from any lateral 
inhomogeneities near the surface so that the f and g values 
along the boundary correspond to those above a uniform 
conductors, Alisio > iff ‘any non-horizontal tdigedcntinulty inter— 
sects the external boundary of the mesh, the distance of this 
intersection spoint sfrom thessurface misit bévlange téenough cso 
that the field values may be relatively small there. In 
the case of a uniform conducting half-space the field values 
are less than about five percent of the surface values ina 
depth of three times the skin depth. Thus, REA is chosen 
EOeb ei a2ius as % a depending on the geometrical distribution of 
the tegions iof “diifernent conductivities iwhere He is the skin 


depth tim the material of the highesteresistivity in athescon— 
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The grid spacings cannot be chosen atbitrcaraiy*® for 
the following reasons: 

(L)eSometrestrictionris) setoforithewadjacent®grid 
Sspacingsyin the sradial ascwellfas in thetangulati direction 
by the fact that finite differences are deed FOr EheCfirst 
derivatives .of ithe fieldsquantities in equations *®(226)P8'This 
implies »that Ar, mustinotedififer fron Ar, very much and 
similarly, Ado, has to be nearly as large as Ad. (See Figure 2). 

(2) The «radial ‘distances Ar, and Ar), may only 
ghichtiiy fdiiiexr fromhthe ehéngths od ftthe tarés RAG, and RA®, 
tee achipridpoimt AWerorcexanp le | \PETthetgrid spacings #4 long 
the sradiimarevgerteralliyemuch*larger thansthée spacings alous 
the arcs, the iteration proceeds separately to a certain 
emtient foraeach elrele, jHowever, ‘the titeration can be per- 
formed correctly even if the ratio of the radial and angular 
grid spacings cannot be chosen properly. In equations (2.8a 
and bp) the “points closest to the point regarded do not have 
tO be used, lr the next neighbouring, point is still in (the 
Same conductivity region, it is possible to take its value 
Grate po even further away from polut Olin Pigure: on 11 
the angular distance Ad does not considerably increase. 

(3) Relatively smaller Ad values have to be used 


mear the surface to obtain More Gorrect He and E.. values 


along the surface (See section 3.3). 
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lt is not possible to obtadnea grid which fillacall 
the above requirements, Jit ds particularly difficult to 
find the proper ratio for the radius of the smallest are to 
the radius “of the mext are, ~ Lf Ar. is Very smalij ali -tne 
Ad values should be relatively large to fill the second 
Fequirement. Also, it is necessary to use a great number of 
radial “erad spacings to Set ee large enough so that the 
titst tTequirement is mot violated, in this work Ar. is 


chosen to be at least three times as large as Ar in this 


Lass 
way the F value at the central point only slightly influences 
the values along the are of smallest radius. The F value at 
this unique point of the polar grid is taken as the average 


of the values at the adjacent points on the left and right- 


hand side after each step of the iteration. 


Sooeeeca cubation of the wSurtiace Valawes 


The surface vaiues of the electric components) are 
Calculated uwsing equations (2..2a)F and ((22.2b) in’ the 
H-polarization case. Substituting finite differences for 


the derivatives we obtain approximate values 
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fhe] Comer, (UNS) 
Ampl (E_) = — 
18 ko , Z 
(re "sti n >) 
(3. 4b) 
oe 
Phase (EL) = 2hetvan G ®? 
ine thie tease ae = 0 all along the surface since H = 50, On 0) 


does not change there, 
In@the E-polartzatton case the magnetic Preld= values 
along the surface are determined from equations (2.1la) and 


(2,16), Using finite differences again, we obtain 


Z Z 
(A£)° + (A_g) 


iL 
Ampl (H,), = —— 
BG (Ar)* 
@io5a) 
-A f 
Phase (Hy) = aretan ( a 
: (Ayf)” + (Aya)” 
Ampl (HD ie) tae a7 = ae ee 
In (2.6 46.n Ad,) 
(345D) 


ad 
Phase (H _) = arctan eee 
ie TO gs 


Sinhd, was taken for Ao, as if the Cartesian co- 


ordinate system were used. 


In the E-polarization case the surface values must 


be corrected by interpolation for the point where the appro- 
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priate arclintersects thevsurface. | A¢ddtnear interpolation is 
used. A region of the mesh near the surface is shown in 


Figure 3, where 


OF =r" cos(A®, ) 
os = ( a t Pi)/2 Goer) 
OF = 1 


The finite differences used in equations (3.5a) and 


(3.5b) sarefobtainednforhthestivvaltecs 


(3.6b) 


and similarly for the ¢ values. 

In the H-polarization no interpolation is necessary. 
E values are determined at point Ce by using the values 
Or & and) pf ab point Py and at the surface. 

The calculation of the surface values from the 
finite differences implies the assumption that the fields 
change linearly downward within the conductor. However, the 


soiution of equation (2.5), in yany. .ayered conductor is) ai athe 


form 
F ¢=ieD exp(-fi nr Sumi 0 pre Dd, exp(vi nr Sin ©O)) G3. 7) 


Where 7 teters te the actual layer. 
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If °therevis no¥verticalediscontinui ty enear ithe 
surface, the upper region of the conducting half-space may 
be treated as a layered conductor. In a polar mesh the 
uppermost*layér is substituted byrtherfirsthsection-of vAd 
nextstoethetsutrrace MethesthicknesseofAthias layer decreases 
as it approaches the centre of the grid from the outer 
regions, At®™thePpecentral* point ity becomes zéronand) thes sub-— 
stitution is not acceptable. Near the extremities of the 
mesh more accurate results are obtained for the surface values 


of He. and E by using the exponential approximation rather 


than the linear one. From equation (3.7) we obtain 
OF : 
Pre cy) Berek Se) ab) (3.8) 


where DS and D, are determined from the field values at the 
surface and along the radius just under the surface at a 


distance of r= sinhd, . The value of n taken here is the 


value of n corresponding to this,section. Also 
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If the values of ae and ane in equations (3).4b) 


and (3.5b) are small, more precise results may be obtained 


using the f and g values at gridpoints that are not those 
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nearest the surface, However, they must be in a region of 
Chemsane SCOnductivity.” in thie case, aii a and E. values are 
calculated from the F values along the radius with slope Ad” 
when Ad” is larger than Ad, -¥ Similarly, the surtace values 
of He can be determined using greater Ar distances in 
equation (3.5a). 

The apparent resistivity values (Cagniard, 1953) 
at the surface are calculated from the values of E and H. 
In the H-polarizgation case the’currents flow perpendicularly 
to the strike of the two-dimensional model. The apparent 


Péesietivety is 


Aa Sas pe PAGS S 
oO. = 7 ~ (—) ‘ (3.104) 


In the E-polarization case the currents flow parallel to the 


strike of the model and the observed resistivity is 


SE NL 72 
0 = i 7 GF ) : (3.205) 


The surface values obtained from equations (3.9a 
and b) differ slightly from those obtained from equations 
(3.4b) and (3.5b) respectively. The only exceptions are the 
small slope models in the E-polarization case (See Fitpure 10). 
Figure 5 shows the amplitude and phase values of H and the 
amplitudes of the apparent resistivity for the slopes 


tan ons On 2 san d= tat a. = o , The curves noted by stars are 


ea siete eR 
abdiay eaettys oid ance, | ls asd 
wh asonaaabb +A rsteeke auten siaaaaatia ad BD 4 


(Geet ,bisineed) 240tev yttubsetess omozaege adit 
‘2 bos & to esnley sds mot? beinivoiss 218. costume afd 38 
yfyalvolbnegrag woli, erasttus ant sexo notiesizetoqg-H eft al 
ebom imnokensmib-ows 443 te yabase saa ge 
er qtivisenaes 


gneisq¢qa sat 


(n0t.) ; (9) ya. ii “9 


eH? 64 [siiletsq wolt etastiss off Sas notissizsLoq=4 edt at . 
ak yiiviseresy boyteado aia bis Lebom ofa to ares, 


eae “o 
ep ates ie 7 


(d0L.&) 


s®.£) enokispps mo3} HaABea Ho @aulsv S562108 of 

anokssups aon} tihivabee eeodt mort yisdgtia 1st2tb ‘ada 7 

ois 976 atiolsgeaxs «ino edt) .vfavisasqaes (AGiE) ba COED 
(8 ssuatt 992) geno gotasshaniog-it, ad ab ehebow aqole tt ne 

SN Pe ee seadqg box shvyitqas eda awole ¢ om 24 7 


sc angele out 2otuttbvidehsss jnadéacs oda 20 


oe wieae yd bazom covsw> ait. © = .9 aay ban $.0 *y 


Figure 5 


29 


HR_ (MKS) 
300,00 


275.00 


250.00 


50.00 -40.00° -30.00 -20.00 ~10.00 0.00 
HORIZ DIST (KM) 


-40.00° 


2% 5¢ 


-45.00° 


-50.00° 


PHASE HR 
-55.00° 


-S0.00 40.00 -30.00 -20.00 10,00 20.00 30.00 40, 00 SO, 00 


~10, 00 0,00 
HORIZ DIST (KM) 


APPRES (OHMM) 
60.00 80.00 100. 00 


40.00 


80.00 40,00 -30,00 -20,00 10,00 20,00 90,00 40,00 0,00 


-10,00 0,00 
HORIZ DIST (KM) 


The surface values of He and 0 “for the sloping 
contact models with tan a, = 0.2 (curves 2) and 
tan a. = © (curves 5) in the E-polarization case. 


The curves marked with stars are calculated from 


equation (3.5b), the curves without stars from 


equation (3.9b)". 
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calculated from equation (3.5b) while those without stars 
from equation (3.9b). The differences between the corres- 
ponding values increase at the extremities of the surface. 
Near the central) point. the «two .cumves, deviatei from cach .ther 
for the smaller slope only. . Curves 2 and 2a show, especially 
for the apparent resistivity, that the exponential approxima- 
tion does not work as well as the linear approximation for 
this, panticular, model,. » Over, the modell rf: thes wer tiical 
contact, however, there is no observable difference between 
curves 5 and oa in) tt hie, ecentrad: wiegiion . omint ttivet H-po Lbarization 
case the E values for all slopes determined from equation 
(3.9a) are acceptable as well as the values from equation 


(3...4b,)+. 


3.4 The Iteration Procedure 


Weert te was, menmtiltomedi wv nt tse citisomsh 2312 fanditr. ito tehe 
calculation of the unknown field values employs the Gauss- 
Seidel iteration method. After the initial values over the 
mesh are set, the iteration procedure continues until all 
the differences between the most recently calculated fo and 
Bo values in equations: (2.1la’ and) b)* or :(2.03a. and b) and 
the corresponding values from the previous step of the 


iteration are less than a given residual € . For each model 
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Chemvialtie®Gf €) istitaken® ast OF 000145 omnis! 1.90.09 percent 
of the normalized f value at the surface on the left-hand 
side of the model. 

The point-by—point iteration proceeds in’ a clock- 
wise direction around the ares of the mesh, neeahnane with 
thesoutermost cirele* and moving inward: *Sintce in’ €ach step 
of this procedure’ the most recent available values of f£ and 
Bo are®used at’? the’ fourSneighbouring points around’ the’ point 
regarded, some degree of asymmetry is shown by the final 
field values even if the conductivity configuration is 
symmetrical’ (See* section’ 4.3)% 

Such asymmetry is investigated for the sloping 


contact models. In the H-polarization case the slopes 


With tany O,. = 0.1 and tan a ©o are studied. The slopes 


ih 5) 


with tan a, = OFZP and" tan a, = © are considered for the 
E-polarization case. Initially’ the» region of higher con- 
ductivity is placed’ on the right-hand side of the model as 
teSshowne in’ Fieuret 10S PAfterward this region is shifted to 
the lLeft-hand®side of the mesh.¥+Thése two configurations 
will be referred’ to as the ‘original" and ‘image’ model 
respectively. The standard clockwise direction for the 


iterations around the arcs of the mesh implies that in the 


two cases the iteration proceeds in the opposite direction 
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with respect to the different positions of the better con- 
ducton,..The tadial direction of the iteration is the same 
for both models, 

The final grid values at the corresponding grid- 
Paints of the two configurations which are images of each 
other deviate from one another for each slope. These 
differences are much larger in the E-polarization case than 
in the H-polarization case since in the latter case the 
iteration proceeds over a mesh whose arcs are half circles 
only (See Figure 4). The maximum differences between the 
tteld values obtained from the: "original' and the [image” 


configurations are observed in the central regions of the 


meshes, 
E-polarization H-polarization 

Maximum differences in f: 

small slope 0.0082 0.0007 

Wiel Toe doe Ont ext 0 ,0090 0.0026 
Maximum differences in g: 

small slope 0.0042 0.0004 

verrucae Contact 0.0034 020003 


The relative differences in the f values near the 


surface are less than 0.04 percent of the f values at the 
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corresponding points in the H-polarization case and less than 
i percent in the E-polarization Case. The relative differ- 
ences in the g walues are much larger. Although in the 
investigated cases they are under 18 percent of the corres- 
ponding g values at the points near the Sey a in other 
cases this value may be exceeded since g values are near zero 
elose to the surface. 

The differences described above are negligible with 
respect to the sSurftace—values—ot E. in the E-polarization case 
and H. in the H-polarization case. However, the surface values 
of the other field components calculated over the 'original' 
models exhibited in Figure 10 deviate from the corresponding 
values determined for the ‘image' models. For the purpose 
of comparison the profiles for the 'image' models are 
reversed in Figures 6 and 7. 

Figure 6 shows the surface values of the electric 
field for the two slopes investigated in the H-polarization 
case. The amplitude and phase values of E. are calculated 
from equation (3.9a). Above the poor conductor and moving 
toward the central point the curves determined from the two 


different configurations deviate from each other for both 


slopes. 
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The surface values of E. for the sloping contact 


models with tan o.0= 0.1 (curves 1) and tan a 


0O 


(curves 5) in the H-polarization case. The curves 


marked with stars are determined from the '‘image' 


model and the’ curves withoutsstats from@the 
‘original’ model as explaineduinwsecti6ns334¢ 


Ce, 


: 8 
ay a 
=i 
8 ° 
; 8 
$+ # 
8 . 2 
;. 8 
gt 9} a 
2 ; Z 
8 o, 
2 8 
A “3 
3 g 
= z 
z 4 
= 8 
Ls : ie sR eI har hi —— 
io toe ‘Cris ethyen the bh oeetries Honest amy OD wen 
8 
a 
8 8 
: 8 
‘| gs 
: 
‘ “. 
8 = 
zg Be 
—8} as 
3 
8 8 
x F 
8 8 
240,00 40,00 S000 
Figure /7 The surface values of He and ES for the sloping 
contact models with tan a, = 0:2 (curves 2) and 
tan a@. = © (curves 5) in the E-polarization case. 


The curves marked with stars are determined from 
the ‘image' model and the curves without stars from 
the 'original' model as explained in section 3.4. 
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Figure 7 exhibits the surface values of the magnetic 
field components for the two slopes examined in the E-polari- 
zation case, The amplitude and phase values of He and -. are 
Calculated rrom® equations (3295) and (3.54) respectively. 

The differences between the profiles fer ehRieen from the 
"original* and the ‘image’ models increase for the He compo- 
nent moving toward=the central” point from both sides for 

both slopes. He values determined from equation (3.5b) would 
show the same features. The as profiles differ only slightly 
from each other above the vertical contact, but there is some 
deviation between the curves from the different models for 
the small slope, especially in the phase values. It is 


apparent that the H, values change less when the direction 


6 


of the iteration is Changed” around” the’ arcs onlys>” since they 


are determined from the field values along the radii. 
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CHAPTER 4 RESULTS FROM THREE TWO-DIMENSIONAL MODELS 


4.1 General Descriptlon of the Results 


The perturbation problem was studied for three 
different two-dimensional models using the described method 
in the cylindrical co-ordinate system. The models. are 
iilustrated in Figures 8, 9 and 10. The first two models, 
the ridge and the dyke, are symmetric while the five different 
slopes in the third model show asymmetric conductivity 
conitigurations. 

Both H and E-polarization cases are investigated 
for each model. The mesh used in the H-polarization 
corresponds to the bottom half of the mesh used in the 
E-polarization for each model. The grid spacings applied 
in the individual cases are presented in Table 1-4 and will 
be discussed later. 


The field values are calculated throughout the 


whole region! of the grid. In seneral 
F = (f + ig) exp (iwt) C41) 
where F = H_ or F = E.. depending on the case being considered. 
Zz 


The amplitude of the field is determined at each gridpoint 
2 
etme mee ennde ce 6 (4,2) 
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Figure 8 The ridge model. 
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Figure 9 The dyke model. 
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Figure 10 The sloping contact models. 
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The amplitude distribution of the magnetic or 
electric field depending on the polarization case is 
presentedirover ‘avlimi ted région tof thermesh wherdathe 
inhomogeneities are located for each model at each frequency 
Investigated, According teteqwa tion Qeleiicin the H= 
polarization case, the contours of equal He represent tthe 
lines of force of the electric field and consequently the 
inves. fof Scurrent “flow “perpendicular oto "the isitrike of the 
two-dimensional model. In the E-polarization case, equation 
(282e) Sshows *that the contours fofveqia l E. correspond to the 
Vines %of “force Sof “the magnetienfield. \cAlsoptatitentmulrei-— 
plying by the conductivities of the corresponding regions of 
the model, these contours represent the distribution of the 
cCuprent “density tilowing ‘paralbleletosthe.strike of the 
SEructure. 

The instantaneous field values observed are 
Galeulated fromtequation ¢C4e1L) 


F = Re[{F] = £ cos(wt) - g sin(wt) R Caos) 
obs 


For each model the field configuration is presented for 
certain epochs of the firet half of the oscillation period 

Tt =217/W over a region where the uniform source field is 
mostly perturbed. During theiyctle the source field: changes 


direction. The configuration of the induced electric currents 
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iS gihie, same at fiilime «t u= s0 ea’s Gat atid= OT / 2 except for the direc- 
tion ,of .the -currents’. » During ithe second ‘half of the “pertiod 
the ,cumments ~ilow «inl the iddrection opposite stot that of tthe 
Pie Sit) Ania iif. 

in tee, H=polarige tition: ic ase thie surface values of 
bhe ~honizonta hveleotric faeldicand) the tappanent rests ivi ty 
from equation, (3,.10a) are presented. The) profile of the 
amplitude of P- Ley SaimisMare ito! thate of E since A is constant 
along sthe surface). ta ths: s uficie Lents ttoivexhibat) onl y\-che 
phase profile of E . In the E-polarization case the surface 
Wranl eS, .0Od st ies erect ree ifiteiid and theshorizontal end® vertical 
magnetic components are determined. Also, the apparent 
resistivitys, from equation’ (3.2 0b)cands thet ratiot of! thes two 
magnetic field components are presented. 

Sdmee, in the E-polarizgation!lcase® the Glectric® field 
and in the H-polarization case the magnetic field is norma- 
lized,’ the-other:£ieldy componentsy calculated fromytheve 
normalized E and H values show relative values. The apparent 
resistivities represent real values. 

The computer programs work in electromagnetic as 
welleas ineMKSeunitsweathe calctiatienstares done ineuks 
units for each model except for the dyke model. 

E and He values.areocalculatediusing the Linear 


as well as the exponential approximation as explained in 
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Seecvlowrts £32, pele evalucado ft dD, and D, in equation (3.8) are 
determined from the F values along the surface and along the 
radius whose vertical distance from the surface is r * sinAd 


Whoreehoni go2es tian 65% 


. Since more correct surface values 
are obtained from the exponential approximation iy the toeter 
region of the mesh and the results of the ‘two ‘approximations 
do not observably differ from each other near the central 
point in most cases considered, the profiles presented in 
this chapter are determined from equation (3-948) or (3.9b) 
depending on *the polarization case. The only exception is 
the sloping contact model in the E-polarization case as 
discussed (in ~wsection 3.3. « 

Each of the profiles shows a gap above the central 
Point Of the mesh, There are two reasons ~10r not retaining 
Che values there. (1), The central point is in a specival 
position. The field values there may be determined using a 
Cartesian co-ordinate system. However, these values usually 
dosnot fit the corresponding profiles since they are calcu-— 
lated from gridpoints whose positions relative to the central 
point ,are.quite different from the relation between the 
points,of the polar mesh.pyi2) The radius of the inuermost 
arc is at least three times as large as the grid spacing 
between this arc and the next one. In this way the substi- 


tution of finite differences in equations (2.6) is only approxi- 
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mate for the points around the circle of smallest radius. 
This may distort the field values over an even larger region. 
painee other authors have published profiles Hin tie 
Cartesian ‘co-ordinate “system; the directions u0of the ‘radial 
andtangular fiieldlcomponents!ares¢orrectedehere so that E. 


and He ape directed to the right and Hy downward jail along 
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4.2 The Ridge Model 
“Wee. General 


A ridge-like model with resistivity 1002m surrounded 
by a material of 108m as shown in Figure 8 is studied. The 
eOpeet the estructuresateitsscenterviseat-d edepth of 23.3 +km. 
The perturbation problem is investigated for two different 
oscillation periods of the source, 10 sec and 100 sec. The 
skin depth in the covering layer is 5.03 km for the shorter 
period and 15.92 km for the longer period. The sets of 
erid spacings used in the mesh are given in Table 1. “Ihe 
fadius of Ehe outermost arc of the mesh is 66.7 km for “the 
LOMsec pertod and 155,35 km for they 100 sec period.) ihe 


calculations are done in MKS units. 
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Grid Spacings for the Ridge Model 
(The radial grid spacings are in kms from the central point 
toward the external boundary, while the angular grid spacings 
aresinvdegrees in a clockwise direction beginning at the 
surface on the right-hand side of the model) 
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4.2.2 H-polarization 


Figure 11 shows the amplitude contours of us tO 
both oscillation periods studied. It can be seen that the 
currents are concentrated above the ridge. Inside the 
lower conductivity structure the current eee spread 
downward, This tendency is better observed for the shorter 
period since the skin depth is less. 

The instantaneous field configurations for the 10 
sec oscillation period are.illustratedsinijFiguréeel3iruAt the 
beginning of the period there are strongly curved negative 
current lines above the sloping sides of the ridge well 
DelowetWes Sit a cena DUttne-—tne Lirst. quarter of the oscillation 
period the currents penetrate deeper into the conductor. The 
current lines near the surface on the left-hand side become 
separated from those on the right-hand side. At the beginning 
of the second quarter-cycle the currents are mainly concen- 
trated above the sides om the!) ridgeyin the better conductor. 
During this quarter—-cycle the; curved “current™lines move doawn— 
Ward and decay and currents of opposite direction increase 
mear the siurtace, 

Higure 14 shows the surtace values of the electric 
field and the apparent resistivity... The profiles feature 
similar changes for the two different oscillation periods. 


However, the increase in the E and p. values toward the 
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Figure ll Ridge model, H-polarization. Amplitude 
contours of Une 


Figure 12 Ridge model, E-polarization. Amplitude 
contours of E 
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Pigure 13 Ridge model, H-polarization. Contours of 
equal H. at equal intervals of one-sixteenth 
of the period. (T = 10 sec) 
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central point is larger for the longer period since a greater 
proportion of the currents penetrates into the region of 
higher resistivity. There is a small dip in. the curves 

above the top of the ridge for both periods. The maxima 
arerlocated=slightly-nearer—the central point, than the. top 
edges of the ridge. “These maxima occur at the places where 
the density of the current lines is greatest. Above the 
middle of the ridge the ES and correspondingly the p- values 
decrease since the higher conductivity layer Gs thicker there. 
The profiles correspond to, those of a,unitonmentedldscuts ide 

a distance of 20 km from the center for the 10 sec period and 


50 ikm for the 100 sec period. 


4.2.3 E-polarization 


Figure 12 shows the amplitude contours of E. Cor 
both oscillation periods studied; It ads apparent for the 
shorter period that the density of the currents flowing 
perpendicularly to the plane considered is much higher above 
the slopes of the ridge than inside it, especially when the 
different conductivities are taken into account. 

The instantaneous field values at eight successive 
epochs during halt of tie cycle are precented In Pigures” 15 
and 16 for the 120 see ogceil lation. period. §At the end cf the 


first quarter of the period the region of maximum current 


' pie en as an 
eevius edd ot qth tinna ab stadt .ydivEvetess 2 
émixem oiT .2holisgq cued ‘yo? ogbbr sda Yo qo aun a 
qo3 sia asda tatoq fetinso sis tsts£et yisdgite haserok 18 
stardw esssiq sid Js tu390 amixsm seonT .sgbia oft re 
adi svodA 229798578 at enti Jnor1ru9 sa Yo ane 
zsulsv 59 gdo vigntbroqesrto> bas 43 eft sg6ts sft Yo —) 
sexent yvedotds et reysl vwiivistoubaos rorgid siz sonke sanatoeb © 

shiejuo bieit mrotinu s to se¢odd of baoqesrzros asiliorg oAT 
bos botesq 952 OL sft to? r93099 add mort md OS to sometatb a - 
-bokrsq 292 DOL ois xox af 66 : 
noisesttalog-% €.8,6 — 
tat <2 fo azuotaos sbuatiqus ad? awote Sl stuglt rh 
aqt <ol dnsikagqe et 31 «betbute ehoixsq aotasiltoso ased = 
gatwolt asastivy sd to vitenab of? asds bobisq t9estosle 

syode szsigind doum et borsbtenos eneigq sia of etyelustbasqreg 
6d4 msdw yviisiosqes ,if sbtent gsit sgbix sd to esqola oft 
sidmuosss o3at qoded avs estitvisoubros Ing79222b- 

“ sviseosouve tigts gs esulsv brsti suoonpinsiant sdT 

af aomugtt nt wtih Scan AC ccd att, Lo 2ae gaizub — 


“nlt 


ous to wil ost 3A. -pouaag got aatttsan 248 or ans tot are 
“auiasres musta to . td botisq “y Po retteup : 1 
re angie ‘o . nt . i i 


50 


Figure 15 


Ridge model, E-polarization. Contours of 
at equal intervals of one-sixteenth 
of the period for the first four epochs. 
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Figure 16 


Ridge model, E-polarization. Contours of 
E at equal intervals of one-sixteenth of 
the period for the second four epochs. 

(T = 10 sec) 
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density lies directly above the ridge. During the second 
quarter of the period this current concentration divides 
into two parts and separate current wedges are formed 

above the slopes of the structure. These current concen- 
trations move downward and decrease in Hare Suge with time 
while the field increases in the opposite direction near the 
surface, 

Figures 17 and 18 exhibit the surface values of 
the fields. The electric field values increase above the 
ridge. The maximum is flat and narrow for the shorter 
period. However, for the longer period there is a dip above 
the top of the structure and the maxima are approximately 
above the edges of the ridge. In this latter case the 
currents penetrate more deeply into the higher resistivity 
region. There are large maxima in the He values on each 
side of the structure for both periods, but these occur 
farther from the center of the structure for the lenzer 
period since the skin depth is Larger. Also, for the 
longer period, some small secondary maxima appear near the 
Central point. Since the ridge is symmetric, He is zero 
above the middle of the structure and changes sign there. 
Furthermore, the amplitudes of H decrease smoothly approach- 
ing the central point and the apparent resistivity prorite’s 


show Similar features to Eo The curves son vthe ratio of the 
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magnetic components are similar to those of pa The apparent 
resistivity profiles correspond to those of a uniform field 
outside a distance of about 30 km from the center for the 
shorter periodt ands 70) kmafory the) longeire perdoda: This implies 


that the ridge affects the field configuration in a larger 


region in the E-polarization case than in the H-polarization. 


4.3 The Dyke Model 
4,3.1 General 


A dyke-like structure with resistivity 4m buried 
in a material of 10082m as shown in Figure 9 is investigated. 
the top of) thet3. 64m, thatk etructdre atsits+center.us dat, 4 
depth of 4 km. The problem is studied for two different 
oscillation periods. For the 10 sec period the skin depth 
Le >..o kmeinrthe poor conductor sand 3/72 km in sthe better 
conductor. For the 40 sec period the skin depths are 31.8 
km and 6.4 km respectively. ~“ The grid’ spacings’ used are 
presented in Table 2. ‘The radius of the outermost are of 


the mescn i6 64.6 km for the Shorter period “and 209790) km tor 


the longer period. “The calculations “are done in electro-— 
; ‘ natant alah 
magnetic unite where the resistivities yaren4 +910 emu and 
13 


10 emu. 
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Grid Spacings for the Dyke Model 
(lire: radral "rid ‘spacings’ are “im kits from *the “central point 
toward the external boundary, while the angular grid spacings 
are in degrees in a clockwise direction beginning at the 
surface on the right-hand side of the model) 
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4.3.2 H=-polarization 


Figure 19 shows the amplitude contours of H for 
both periods studied. The effect of the dyke is greater for 
the Pshorter-period-—since™the skin depth is less and tHe 
current is denseryin the region jof the structuce.—-—lven for 
the jlonger period the Concentration..of the current -lanes by 
the (dyke whose conductivity is higher—-than-thet—of—-the 
surrounding material is observed. 

The contours of equal H, at different epochs of 
Ene lO} sec joscillation period are exhibited in bvicure 21> 
the current Lines are stromegly retracted at the boundaries 
eofethe structure, During the first quarter of the period 
positive currents décrease and finally current vortices 
form .in—-therbet ter conducting region... During the second 
quarter negative currents increase near the dyke. 

Tressurtace=values of E and Q. trish opunes 22 show 
the asymmetry which was discussed in section 3.4. This is 
observed particularly for the lomger-oscillation period: 
There is a-slight difference between the two periods injthe 
apparent resistivity profiles. Outside a distance of about 
20 km from the center the effect of the dyke is not observed. 
E and p- values show some increase above the structure for 
the longer period. Also, minima are evident just above the 


dyke. The distance between these minima is’ less than the 
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Figure 19 Dyke model, H-polarization. Amplitude 
contours of H. 


Figure 20 Dyke model, E-polarization. <Amptitude 
contours of Bb. 
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Dyke model, H-polarization. Contours of 
equal H_ at equal intervals of one-sixteenth 


of the period. (T = 10 sec) 
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width of the top of the structure. “Since the vertical size 
of the dyke exceeds the skin depth for the shorter period, 
the currents penetrate only slightly into the region of 


higher resistivity below the structure. 


aoe s t"poLari 2? atd.on 


Bigure’ 20 exhibits the tmplitdde contours of E. 
for both periods investigated. For the shorter period the 
current density is higher above and inside the buried dyke 
than at a corresponding depth in the surrounding poor con- 
ductor farther away from the structure. For the longer 
period the contours do not deviate greatly from a uniform 
coniriguration. oince the conductivity is 25 primes higher 
inside the structure than in thelLsurroundine material,-iLt 
is- apparent that the current density increases. inside the 
dyke, 

The contours of equal E. at different epochs of 
the 10 see oscilDation period aré shown in Figure ~23 > The 
eurrents (Cilowinge parallel to therstrikewol the dyke sare 
concentrated inside the structure during the whole cycle. 
The contours, which correspond to the lines of force of the 
magnetic field} are strongly curved near the structure, “At 
the beginning of the £itst and second quatter of the 


oscillation period current wedges are observed around the 


better conducting region. 


oshe Inokitev sits somte vanersuase 903 To goo ada ooo 
vsotzsy Yetxose oly 10? Haqsb athe oda abosoxs aly. fin 
Yo nokgsy sd#3 otmt ylidatiea vino statssa9¢ earners | 
.srudsprte sit woled lacie nl raat 


“ nol seotaloed tek 
= ot 


~ 


EF 


a to etucinos sbustiqms afi exytdidxe ie stygit 
oi? botyeq tesv0e oft rol .beatkgtiesvat ebotieq Hod bial 
si¢b bsixrud oft sbtent bas svode yetigtd ef yolensb snetTy. Y 


-noo toog gatbnuorrue od3 ot daqgesh gotbnogestxo5 6 38 asdd 


segnol sd3 10% .studovrte eft? mort yews t9d348% 203 9ub- *t 
mrotinu & mort yiteery stetveb ton ob savotnos oft bofirseg : 


sergitd asmty @& et ystvidosbnos of4 sonte ,tokistugttmes 

+i ,atxrstem gobbavottue 93 nt aeit sautoutie sid? sbteot 
4dy sbbant esaasiont yttensb tasviv> of? tad3 Jesteqqs ak . 
to efoogs Jnorstith % <i Peups to etvotnos sdT 


eat .€S stugtt nat awotle sie botrsg mottalitoao 298 or ods a 


, 
ots givb sf4 to satiate odd 03 Ielleisq agmiwolt evadrsea “ 
-9foyo siortw sd4 gotz0b stutoutse sit sbtent borsaimesned a 7 
otis Yo soyvoi Yo sent sd3 oF Broqeextos dotdw ,ssvotmes sat 
3A «8. 8IUIouxIe soo tHBOn beviriwo yigaorie sts «bier? ‘jen s 
sHd3 to 'xs3teup booose bes sett? via to gutantged odd 


eit boyows bevisado exe esgbew iasaivo bokasq | 


62 


Figure 23 


Dyke model.’ E=-polarivation.= * Contoure ‘or 
equal E_ at equal intervals of one-eighth 
of the petiod:, 
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Figures 24 and 25 illustrate the surface values 
of the field components. Less asymmetry than in the H- 
polarization is observed in the Ho. and correspondingly in 
the p. values. The apparent resistivity values differ from 
each other=for the two~periods and this difference is much 
greater than in the H-polarization case. The surface values 
outside a distance of about 40-50 km from the center of the 
Structure nearly "correspond to those of awudifrenm Eleld for 
the 40 sec period... This distance is less for the 10! sea 
petiod. The resistivity above the dyke is greater for the 
longer period than for the shorter since the currents pene- 
trate into the higher resistivity region under the structure. 
The change in the phase values of He and p. LS Opposite. tor 
the two different periods. Also H shows a greater and wider 
maximum for the longer period. E. values decrease above the 
structure. The minimum.valuvue, is.less for the shorter period 
since the thickness of the dyke exceeds the skin depth in 
the better conductor. For the 40 sec period there is a small 
local maxinum ti the E. Values above the top of the structure. 


This, shows thes effectsof the, thicker dower conductivity region 


above the center of the dyke. There is a characteristic pair 
CLomaxama dm thie ie values. The distance between them is 
greater than the width of the structure. The maximum value, 


which occurs approximately at the same position for both 
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petiods, is greater for the, donger period. These maxima are 
observed in the Hy fH values as well, but the maximum value 
for the shorter petiod becomes larger than that of the longer 


period. Another pair of maxima in the H, values, which is 


o 

present just above the dyke and only for the longer period, 

ts much less dominant in’ the protule of the tacio of the iLwo 
\ 

components. As with the ridge model, the vertical magnetic 


component is zero above the central point and changes sign 


there. 


4.4 The Sloping Contact Model 
4.4.) “General 


The five- different sloping contacts studied in 
the present work are shown in Figure 10. The resistivity 
Of the Material is 10%im above the sloping discontinuity 
and 1002m throughout the remainder of the conducting region. 
tMesoscl Llatitom perlodmoL the alternating ssource, fica 
40 see in each case. This period gives a skin depth of 
LO.O0) km fou the better conductor and )31..6 km) ton the high 
resistivity material. Ihe two difierent ‘sets of erid 
spacings used in the radial direction=are presented an 
fap le 3. “[he maximum radius of \Ehe mesh ic) 209.5) kmein 
set A and 03.2 km an set B. Ihe foun sets of angular 
Brid spacings which are applied for the ditterent siopes 


are given in. Table 4.) “The calculattons are (done inyMKs> unate. 
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Wailea. 8 


The Radial Grid Spacings for the Sloping Contacts 


(Ay walues are in kms from the central poant Coward 
the external boundary) 


: E ot det 
may: 
eYosdnoD aentgole sda tot egaiosgqe bryd Intbsa ont ft 


9 ami at ste esbisv tA) - 4 


bspwoo tniog fnxtaso9 ons mor 
(visboauod Isatsdxs sit 


Ef 20.0 GRAVE pSal wed aby el, ot. VB. 


Ooh 68 LuN [SB (Ee BLE eRub | EQk Che es 
a 


(Ad 


at 


Set 


Set 


Set 


Set 


68 


Table 4 


The Angular.Grid, Spacingspfiorathe,Slopingr Contacts 


values are in degrees in a clockwise direction beginning 
the surface on the right-hand side of the model) 
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4.4.2 H-polarization 


Inthe H-polarizationti¢ase only four of the stloping 
models are studied. The following sets of grid spacings are 


used for the different slopes: 


B-b for 4, = ee 
O 
B-c for oy =e 
B-c fOr oh = oe 
fe) 
B-c fOr 0. = 90 


Figure 26 shows the amplitude contours of He Lor 
each slope. Tt ads well illustrated! that the electric currents 
tend: to, ilow in thei region of hileher conductivaty. In the 
vertical contact case the current flow depends on the differ- 
ent@iskin depths! in the two regions: | The curpene jJines ‘are 
strongly refracted at the boundary between the different 
GConducthivitty negions. This refirnaction#of the current lines 
is caused by slowly varying surface charges at the interface 
as Kiscuseed\by Jones and, PricexGl9/7iay 19V71ic) \and Price. (1973). 
Lt is, the surface charges builtijup at yany interrace, inside 
the conductor which) spread the lournent wlines vat) thes discon— 
ftinuitiesiin the other two modelstas well... Sinee there is 
more than one discontinuity in the ridge and dyke models, 
the picture is more complicated. In general current lines are 


refracted at all interfaces near lateral inhomogeneities. 
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Ihe field configuration at different epochs of the 
oscillation period for the 45° slope are exhibited in Figure 
27.\ At the-beginning of -the period the distortion ofthe 
current lines, which aréjparallel to the-surfate at great 
distances from the discontinuity, increases with depth near 
the discontinuity, especially in the better cenductor. After 
a quatcrer=period ‘current loops are formed near thegsurtace) in 
the, higher=conductivity region. | Durine the second-quarter 
of the period these current vortices move downward and decay. 
For the smaller slopes the current lines are less curved than 
in the case discussed above. For the vertical contact there 
is an even stronger tendency to form current vortices in the 
better conductorxturine thesescil lation. 

The amplitudes of the surface values in Figure 28 
show a discontinuity at the central point. Near the interface 
on the higher resistivity side E and p- are much greater than 
their-asynptoticoraliiesem Apart from tie vertrealecuntacc., 
whieh Shows-a slight winerease jonly,.—-the-values near phe inter 
face increase with increasing Q@ and the maximum value appears 
Beng Oy = 45 , The surface values correspond to~thosée of a 
uniform field outside a distance of about 30-40 km from the 
center depending on the slope of the discontinuity on the 
Other side. On the higher conductivity side the amplitudes 
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distance from the center where the discontinuity intersects 
the surface. Also, these values approach their asymptotic 
values at a distance where the discontinuity is at a depth 
approximately equal to the skin depth there. The phase 
profile reaches ae which corresponds to a uniform field 
configuration, farther away from the center. For the 45° 
slope the amplitudes near the discontinuity are smaller than 
the asymptotic values and show a minimum at a distance of 
avout 4ekmefrom the ,central point,, The curves correspond to 
thOSe~ Ota Upnittiorm £1eld outside a,ydistancé .of, about.,~25 km 
for the amplitudes and 35: km for the phases, In the vertical 
contact case the amplitude drops to a very small value and 
the asymptotic/value isreached at.a distance of about 15 km 
while the asymptotic value in the phase is reached at about 


25 km from the interface. 


494.3 E-polarization 


In the E-polarization case different combinations 
of the sets of radial and angular erid spacings, are used so 
that the effect of changing the mesh may be studied. The 
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Figure 29 shows the amplitude contours for the 
electric field for the slopes that have been discussed in 
the H-polarization case as well. The change in the field 
configuration depending on a shows the effect of the region 
of higher conductivity. The currents eraeize parallel to 
the strike of the slope are highly concentrated near the 
Ssuriace Inethe better conducting) region- 

The contours of equal E. at different epochs for 
orée=taif ofthe "cycle are exhibited™-in- Ficure=30 forthe 
45° slope. At the beginning of this half of the cycle the 
direction of the magnetic field presented by the contours 
is-—aimost lorizontal except near ttherdiscontinutty._Attec 
One—-erpitn of tie pertod-—cehesadirectilonus-of 7phe=teldecian ves 
and—is—mainly vertical..-Also, in tthis—time ititerval—the 
field magnitude has decreased. During the second quarter 
of the period the magnetic lines of force form loops in the 
better conductor. This shows that the corresponding positive 
currents are concentrated there. There are increasing 
negative currents near the surface. 

The sUrhacée vValues¥calculated” usine jset A for the 
radial grid spacings are presented in Figures 31 and 32 while 
Figures 33 and 34 show the surface values determined applying 
set B. There is some difference between the corresponding 


profiles in the twotséts*of figures#for thestwolsmaliaslopes 
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near the central point. The calculated surface values depend 
on the size of the mesh to a certain extent, patttctitar iy eit 
the maximum arc of the grid is not large enough and there is 
a considerable deviation from the uniform field near the 
surface at the extremities of the model as discussed in 
Secerone3s «2, In the case of the smallest slope the depth 
cebehesdiscontintity “is slightly “eréater than. the ®ekin-depth 
for the better conductor at the external boundary of the 

mesh and the boundary conditions are not satisfied perfectly. 
Also, the iteration cannot be done in an optimum manner inside 
and adjacent to the better conductor. The values of Bo> Bos 
B, and By in equations (2.8a and b) are not of the same order 
of magnitude sinee the distances from point '0' to points "2' 
and '4' are much greaterethan the “distarvees*te points “and 
V3 i€seeehi gare 22): 

Differences in the electric field values are 
negligible except for the greater dips in the middle of the 
curves as calcwlations are made nearer “the ‘central point. 
Differences are evident in the horizontal magnetic field values 
and correspondingly in the apparent resistivity values in the 
@entral region. siThe vertical ’contdet profiles Yarevequivalent 
im Cele stwo sfigwures'. 


The He. and H, amplitude curves both show maxima on 


> 


the better conductivity side. These maxima are shifted to- 
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ward the higher resistivity region as a0 increases. The 
maxima in the horizontal magnetic field values appear farther 
away from the center than they do in the vertical component. 
As a result there are maxima in the Hy fe. values which are 
narrower and shifted toward the center when compared with the 
€urves.o£ He and Hes The maximum for the 26.7° slope lies 
between the maxima of the profiles for the 11.7° and 45° 

slope for both magnetic components and correspondingly for 
Ehevrattosasawell.y The ne component has another even greater 
Haximumvon the higher, ,resistivity.side.. This is quite near 
theypoint where thesdiscontinuity,intersects.the-surface. 
Ltieshifts slightly toward, the better conductor,as, @ increases. 
In the vertical contact case the two maxima coincide at the 
centralepoint, »Sincey the horizontal,component, of the maguezic 
field does not show considerable changes above the lower 
conductivity region, the maximum of the vertical component 


dominates in the ratio H /H . 


0 
The apparent resistivity values show a gradual 
change between the two asymptotic values. The lower resis- 
tivity value is more slowly approached in the case of the 
smaller slopes since the depth of the higher COndUCEIVLEY 


material is less than the skin depth even at large distances 


from the central point. 
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ACcCOvdinegigtahtne surface values, in particular the 
Qe and ee curves since H does not change much, the field 
configuration can be considered as uniform outside a distance 
of about 40-50 km from the center on the high resistivity 
Side. “This distance is larger for the vertical contact than 
for the small slopes. In the H-polarization case this dis- 
tance is the smallest for the 90° slope. Above the slope the 
uniform field configuration appears at a distance of about 30 
km Eor the vertical contact. The amplitude of the apparent 
resistivity reaches its asymptotic value at an even smaller 
distance. This distance increases with decreasing O and for 
the smallest slope it “exceeds the radius of the outermost arc 


Of the mesh,  FRe-discontinuity ‘atfects the freld configuration 


near the surface at the external boundary as discussed above. 
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CHAPTER 5 CONCLUSIONS 


The perturbations of a uniform alternating electro— 
magnetic field by three different two-dimensional conductivity 
inhomogeneities has been studied in this thesis. The cal-— 
culations have been made using a cylindrical co-ordinate 
system which is more appropriate than the Cartesian co-ordinate 
system for the models considered, The field configurations 
near the non-uniform region of the conductive area have been 
given as well as the change of each field component along 
the surface. 

Several sloping contacts were examined by other 
authors using different methods. The present study applies 
a more acceptable boundary sande vonlet han that of Geyer (1972). 
Also, the present work gives a better approximation to the 
geometry of the slope than that used by Jones and Price 
C90 7a) 5 Furthermore, more details about the surface values 
are given than those in the paper by Reddy and Rankin (1973). 
Also, greater slopes than Blake (1970)) uséd are examined, 

Some horizontal inhomogeneities were investigated 
by Jones and Price (1971b) changing the depth and vertical 
size of the structure. However, there was only one oscillation 
period studied. Also, Rankin (1962) examined the electro- 


magnetic effect on a dike just below the surface in the 


H-polarization case only. 
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some limitations’ of thes finite’ difference™method, 
particularly considering the choice of mesh (see section 
3.2) may be examined in more detail in future. The asymmetry 
effect of the iteration procedure (see section 3.4) must be 
further investigated as well. Also, better approximations 
ror the first derivatives of the field values in equations 
(2.6a and b) should be considered. The approximation used 
here assumes a linear change of the first derivatives which is 
inaccurate near the inhomogeneities even if the adjacent grid 
Speacings, do mot difier Erom each other. At the extremities 
of the grid neary~the surface |théerfiteld values *change more %in 
the relatively large angular grid spacings than is given by 
the linear approximation because the electromagnetic field 
decreases exponentially with depth inside a uniform conductor, 

The substitution of finite differences may cause a 
slight change stnmathagtteldevalues: insthée E-polbarlzation *case 
just inside the external boundaries along the surtiace where 
the uniformilfteld cont Peuratiion Cis’ not distorted *by “the 
inhomogeneities. This effect is particularly observed in 
Figure 17 where far tiron. (the ridgectietaurrace valies of "tive 
electric field are less than 1 which is the normalized E. 
value at the external boundary. For the sloping contact the 
electric field increases over the normalized surface value on 


the high resistivity side as illustrated in Figures 32 and 34. 
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